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ABSTRACT 

The interaction between refraction from a distribution of inhomogeneous plasma and 
gravitational lensing introduces novel effects to the paths of light rays passing by a 
massive object. The plasma contributes additional terms to the equations of motion, 
and the resulting ray trajectories are frequency-dependent. Lensing phenomena and 
circular orbits are investigated for plasma density distributions A oc 1 /r^ with h > 0 in 
the Schwarzschild space-time. For rays passing by the mass near the plasma frequency 
refractive effects can dominate, effectively turning the gravitational lens into a mirror. 
We obtain the turning points, circular orbit radii, and angular momentum for general 
h. Previous results have shown that light rays behave like massive particles with an 
effective mass given by the plasma frequency for a constant density h = 0. We study 
the behaviour for general h and show that when h = 2 the plasma term acts like an 
additional contribution to the angular momentum of the passing ray. When h = 3 the 
potential and radii of circular orbits are analogous to those found in studies of massless 
scalar fields on the Schwarzschild background. As a physically motivated example we 
study the pulse profiles of a compact object with antipodal hotspots sheathed in a 
dense plasma, which shows dramatic frequency-dependent shifts from the behaviour 
in vacuum. Finally, we consider the potential observability and applications of such 
frequency-dependent plasma effects in general relativity for several types of neutron 
star. 
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1 INTRODUCTION 


It has been suggested that observations of 
black holes may be affected by general relativ¬ 
ity (GR) and t he frequency-dependent proper¬ 
ties of plasma (iBisnovatvi-Kogan fc Tsupkd l2009l : 


iMorozova. Ahmedov fc Tursunov ~ |2013| ) and that these 
effects may be relevant for black holes i n X-ray binary 
systems llBisnovatvi-Kogan fc Tsupkoil2010l) . Lensing effects 
modified by pla sma surrounding g alaxy scale lenses were also 
investigated bv lEr fc Mad ll2014 1. requiring observations at 
MHz frequencies. In this work we will study ray-tracing 
under the influence of both gr avitational lensing and refrac¬ 
tive p lasma effects following iTsupko fc Bisnovatvi-KoganI 
(l2013r) . and examine frequency-dependent orbits for light 
ray s. Deflection by both gr a vity a nd plasma was discussed 
by iMuhleman fc JohnstonI (Il966li for the case of radio 
waves passing by the Sun, using the product of the grav¬ 
itational and plasma refractive indices to derive a weak 
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field approximation. IMuhleman. Ekers fc FomalontI (llQTOli 
measured the deflection of the radio loud quasar 3C273 
as it was occulted by the Sun and obtained an estimate 
of the electron density in the solar corona. In this work, 
we make use of the geometric optics approach developed 
by ISvng3 lll960l) . which has applicat io n to a diver se range 
of phenomena llDressel et al.l l2009l : iLine 3 l2009li . Fully 
covariant radiative transfe r through a refractive medium 
was treated in depth by Broderick fc BlandfordI ll2003li , 
and iGedalin fc Melrosd ll2001^ discuss dispersion relations 
applicable to ray propagation in relativistic plasmas near 
black holes and rapidly rotating neutron stars (NSs). 
However, even the most basic dispersion effects and simple 
distributions of plasma near relativistic objects provide 
novel examples that differ from the vacuum behaviour of 
rays in the Schwarzschild space-time in interesting ways. 


We review the theory and background of relativistic 
light bending including plasma effects in Section In Sec¬ 
tion O we study the trajectories of light rays through the 
plasma and consider waves near the plasma frequency with 
density varying as an inverse power law with index h. When 
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h — 0, the plasma frequency is constant and light rays be¬ 
have like massive particles with m ass given by the plasma 
frequency llKulsrud fc LoebI Il992h . When h = 2, we hnd 
that the plasma term acts like an additive contribution to 
the angular momentum. The h — 3 case has similarities with 
the behaviour of massless scalar fields on the Schwarzschild 
space-time. Section[3T]includes a study of circular orbits for 
a general power-law index. We apply the analytical results to 
the pulse proHles of a compact object sheathed in plasma in 
Section|4l and review relevant properties of compact objects 
in Section [5] 


2 THEORY 


In this work we use natural units, with G = c = 1, unless 
otherwise stated. As we proceed through the derivations of 
the main equations, we keep the h factors, but the numer¬ 
ical calculations are p erformed with h = 1. We follow the 
metho d and notation of[S vngel lIlQGOlj , and review the deriva¬ 
tion of lTsupko fc Bisnovatvi-Kogan (l2013l l to find the gener¬ 
alized trajectories of photons in the presence of an isotropic, 
inhomogeneous pl asma distribu tion, though more general 
possibilities exist (IPerlickI l2000h . Let us assume spherical 
symmetry and adopt the line element 


ds^ =—A{r)dR + A{r) ^dR + RdO^ + R siid{9)d(R (1) 


with 

Mr) = 1 - ^ ( 2 ) 

and Tg = 2M for the Schwarzschild metric. We write the 
coordinate indices as {t,r,6,(j>), with the index a to sum 
over the spatial coordinates and i,j,k to run over all four. 
Without loss of generality we consider orbits in the equato¬ 
rial plane and set 9 = 7r/2. Now, suppose a distribution of 
plasma surrounds the mass M, which is located at the origin. 
The plasma has an index of refraction n = n{x°‘,uj) where 
the photon frequency is uj, with effective energy Tiui = —piV^ 
with respect to the velocity of the plasma medium Y*. The 
linear momentum of the photon is p* = TiK'', where A* are 
the elements of the wave four-vector. The medium equation 
dlx ngelllQGOll gives the relation between the components of 
the photon four-momentum and the refractive index of the 
medium. 


n 


2 


= 1-f 


(pMM' 


(3) 


In the vacuum case n = 1 and the standard condition for 
null geodesics pip* = 0 is recovered. We describe photon 
paths through the plasma in the Schwarzschild space-time 
by the Hamiltonian in the geometric optics limit. Using the 
variational principle and the medium equation we have 


H{x\pi) 


1 

2 


g'^MiPj + {n 


1) [PiM)" 


= 0 . 


(4) 


Following iTsupko fc Bisnovatvi-KoganI (|2013h . we assume 
the plasma is a static, inhomogeneous medium such that 
Y* = = 0. This static condition simplifies the ef¬ 

fective photon energy Tiuj = —Pt%/—With this condition 
there is no time or </> dependence of the Hamiltonian, so pt 
and are constant, and we restrict pg = 0 so that the orbit 


remains in the equatorial plane by symmetry. The time-like 
component pt is related to the vacuum photon energy 


Pt — klUJoo 


(5) 


and for a static medium this gives 
Pt\/-g** = -pM-gtt = -Tux{x°‘) 
with the effective redshift relationship 
u{r) = 


Ai/2 


we have 

t 

^ = — 


( 6 ) 

(7) 

( 8 ) 


The radial dependence of the photon frequency is due to the 
gravitational redshift supplied by the space-time metric, and 
we denote the photon frequency far from the origin r oo as 
u)oo- For the sake of notational simplicity we will suppress the 
radial dependence when referring to the frequency variables. 

In order to make analytical progress at this point we 
will introduce a specific form for the plasma frequency. Let 
us assume that the refractive index is of the general form 


2 


n 



(9) 


where W = W{x°‘,uj). With equations (O and (|8}, the 
Hamiltonian in equation 0 simplihes such that 


H{x\pi) = i {g'MiPj + fdW^) = 0. 


( 10 ) 


The paths of light rays are then described in terms of the 
affine parameter A along the trajectory by 


dA 


OH 

dpi 


^7 

= g P] 


( 11 ) 


and 

dpi dH 

dA dx^ 


-\glhjPk - 


( 12 ) 


These equations give the positions and momenta of photons 
through the space-time under the influence of gravity and 
the optical medium. The r and (j) components from equation 
im can be combined to eliminate the affine parameter A, 
resulting in 


d<j) _ P±_}__ 
dr pr r^A{r) 


(13) 


The denominator can be simplihed by using equation uni 
to give Apr in terms of the constants pt and p,/,, such that 


Apr 



(14) 


Here we take the positive solution corresponding to the case 
when both the (f> and r coordinates increase, and the negative 
solution with decreasing radial coordinate. To proceed, we 
write the constant p^ in terms of the impact parameter b, 
such that 


P(j) — fiiJccb. 


Then we have 
# ^ 

[l^ (1 - 


A11/2' 

72-J 


(15) 

(16) 
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Gravitational tensing in plasma 3 


As expected, for W ^ 0, equation m reduces to the 
vacuum case. This equation gives the paths of photons in 
the equatorial plane in the general case of strong deflec¬ 
tion which is appropriate for orbits near the surface of the 
star. The refractive index of the plasma in this e xpression 
intro duces frequency dependence to null geodesics JNoonarJ 
I 1982 I I. With the substitution u = M/r, outgoing rays have 


dij) 

du 


M2n2{u) 

52 


1 


(1 — 2u)u^ 


1/2 ’ 


(17) 


where the potential dependence of n on u is explicitly stated. 
Now we consider an observer a great distance away at r = 
D —^ 00 , such that the fiducial direction from the centre of 
the star to the observer corresponds to the impact parameter 
b — 0. Since we are interested in describing rays which leave 
the surface of the star {upt = M/R) and reach the observer 
(ud = 0 ), we define 


A = 




du 


1/2 ’ 


(1 — 2u)u^ 


(18) 


with the additional simplihcation x = h/M, and taking the 
sign into account through the integration limits. This ex¬ 
pression gives us the maximum viewing angle on the stel¬ 
lar surface that can be seen by the observer. In flat space 
Amax = 7r/2 always, which means exactly half the surface of 
the star is visible. This is not the case in GR since gravita¬ 
tional lensing acts to redirect light rays such that a distant 
observer sees a larger apparent area and more of the stel¬ 
lar surface than would be possible in flat space. In fact, in 
extreme cases the observer can see portions of the surface 
multiple times (for Amax > rr). 

It is also useful to describe the impact parameters of the 
rays which connect the surface of the star with the observer. 
Let us define 5 as the angle between the radial and angular 
components of the photon momentum. At the surface of the 
star r = R, this is the angle with respect to the surface 
normal that the ray is emitted 


tan 5 = 


\/p^ 


y/P^Pr ' 

Using equation (HI leads to 
Rn{R) 


A{Ry/^ 


sin S 


(19) 


( 20 ) 


where the maximum impact parameter fcmax occurs for 
5 — -E12. This formula reduces to the standard expres- 
sion for the emission angle in vacuum when n = 1 

jPechenick. Ftaclas fc Coherilll98'^ : lBeloborodovl[2002h . The 
presence of the refractive index within the impact parame¬ 
ter expression predicts that the portion of the stellar surface 
visible to a distant observer shrinks as n —>■ 0 . 

Specihc numerical examples require a choice of 
W{x°',uj). The simplest option is to pick 


= ujI 


AEC^Nir) 

m 


( 21 ) 


where uie{r) is the plasma frequency, e and m the electron 
charge and mass respectively, and N{r) the plasma num¬ 
ber density. Together with equations (0 and 0 , equation 
(HB) gives the dispersion relation for rays in a cold, non- 
magnetized plasma. We will make use of this form in the 


following sections, but we stress that our results to this 
point are general for n = n{x°‘,uj). Propagation through 
this medium requires that the photon has ta(r) > cae(r) ev¬ 
erywhere. In the case when uo ^ uje ox the plasma density 
vanishes N{r) = 0, we recover the vacuum case n = 1. The 
condition for rays to propagate from the surface to infinity 
is 


CUoo > LO,{R)A{Rf/^. 


( 22 ) 


Consider the situation illustrated in Fig.[T] In this hgure we 
show the behaviour of light rays with a range of asymptotic 
frequencies ujao- We use the surface value u)e{R)/uj{R), with 
equations 0 and the propagation condition to label 
the frequency ratio. Tracing rays from the stellar surface 
numerically requires a specihc plasma density distribution, 
so we choose a physically motiv ated example based on th e 
Goldreich-Julian (GJ) density llGoldreich fc Julian! I1969I L 
Normal pulsar magnetospheres contain plasma with densi¬ 
ties that are sustained by GJ currents, transferring charge 
carriers from the surface of the NS to the magnetosphere. 
These currents therefore provide a lower limit for the den¬ 
sity of charge carriers in a standard magnetosphere. Using 
the plasma frequency (equation [ 2 TJ, let us consider a radial 
power-law density 

N{r) = (23) 

where h> 0, such that 



(24) 


Since the GJ density depends on the strength of the polar 
magnetic held we use h = 3 giving a 1/r® form. Let us 
consider a highly relativistic star with i?/rg = 1.60, and 
assume a plasma frequency cue with fe = 1 as an arbitrary 
constant. We plot the rays with positive impact parameter 
ranging from 0 to 6 max as thick dotted lines. All frequencies 
satisfy the condition for propagation given in equation (I22II . 

Classical turning points exist for the ray paths where 
dr/d(f) = 0. For all rays with impact parameter b < fcmax 
these turning points are within the spherical surface of the 
star and have no effect on the ray paths. For rays with the 
maximum impact parameter the turning point is exactly on 
the surface of the star, indicating a ray that just grazes the 
stellar surface. The rays which have larger impact parame¬ 
ters than 6 max do not connect the stellar surface to a distant 
observer. 


3 NULL ORBITS WITHIN REFRACTIVE 
MEDIA 

In this section, we are interested in the gravitational lens¬ 
ing effect on rays passing by the star. These rays approach 
from infinity and have impact parameters b > 6 max. The 
turning points for these trajectories, Vp, are found when the 
denominator of equation (HU vanishes: 



Using a power-law density distribution gives rise to a poly¬ 
nomial with order that depends on the power-law exponent 
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Figure 1. Examples of ray-tracing from the surface of a compact 
object for a sample of photon frequencies, denoted by surface 
values Ue{R)/Lo(R). The vacuum solution is shown on all panels 
as a set of thin dashed lines. Impact parameters between 0 and 
^max Q-re used to illustrate the shrinking of the visible surface for 
rays that approach the plasma frequency. This example makes 
use of a highly relativistic star with R/vg = 1.60, and a plasma 
frequency of the form = 1/r^. 


CD (R)/co(R) = 0.9 





+ 2M-|- = 0. (26) 

ojio 

The solutions of this equation signify the boundary between 
regions in which A is real and complex valued. A trajectory 
that reaches a turning point has vanishing radial momen¬ 
tum and subsequently reverses direction. When the effect 
of refraction is significant, even light rays near the line of 
sight with small impact parameters can be strongly affected 
by the interaction with the medium and turned away. In 
this case the plasma behaves like a scattering medium, and 
it is only far from the origin that the vacuum-type lensing 
behaviour is recovered. 

We define the energy and effective potential in analogy 
with the vacuum case by associating the with the angular 
momentum of a photon, and the time component pt with the 
vacuum energy of a photon. We call the radial rate of change 
with respect to the curve parameter f, and equations (El) 
and m can be used to write 



The effective potential of the Schwarzschild space-time with 
a plasma distribution is 




0.2 

0.15 


'lu 0.1 
0.05 

0 

0 5 10 15 0 5 10 15 

r r 

Figure 2. Ray—tracing and effective potential, h = 1 case: an 
example of ray-tracing in the case of u}^ = k/r for rays that 
have u}e{R)/u){R) = 0.99, for R = 3.2 and M = 1. The potential 
acting on the dashed ray is shown in the lower-left panel and the 
potential acting on the dotted ray in the lower right. The effective 
potential is the thick line, the vacuum and plasma contributions 
are the thin and dashed lines respectively. 




where I4ft = V + Vp, with V and Vp the vacuum and plasma 
contributions, respectively. In the absence of plasma I/p = 0 
(fe = 0) the potential reduces to the vacuum case for null 
geodesics. 

For the h = 0 case oie is constant, and photons be¬ 
have like massive particles with effective mass equal to 
the plasma frequency , a scenario that was explore d by 
KuIsrud^^jO^ (Il99d~). Broderick fc BlandfordI (l2003l ~) and 
Tsupko fc Bisnovatvi-KoganI (l2013h . These photon paths are 
identical to the well-studied time-like paths in the vacuum 
Schwarzschild case. However, it is far more interesting to 
explore the behaviour for other, less trivial plasma density 
distributions. The h = 1 case is plotted in Fig. [ 5 ] for the 
trajectories of 10 rays past the mass M = 1 with radius 
R = 3.2, and shows the corresponding effective potential for 
two of these rays. We plot the square of the photon effective 
energy — pi, the effective potential I4ff, and the contri¬ 
butions from the vacuum and plasma parts V and Vp. The 
effect of the plasma dominates over gravity for rays near 
the plasma frequency, and this scattering behaviour is more 
pronounced at lower h. For the h — 3 case, shown in Fig. [S] 
the plasma density drops off quickly enough that we recover 
behaviour similar to lensing in vacuum for rays further from 
the compact object even at low ratios of photon to plasma 
frequency. In the figures we arbitrarily set fe = 1 for these 
numerical examples. 

The h = 2 case has the plasma term act like an ad¬ 
ditional contribution to the angular momentum. To make 
this connection explicit, consider the turning points from 
equation (I26II with the angular momentum L — ojoob: 




(29) 
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r r 


Figure 3. Ray—tracing and effective potential, h = 3 case: ray¬ 
tracing with uJq = k/r^ for rays near the plasma frequency 
LJe{R)/u{R) = 0.99, and other parameters as in Fig. [2] Note that 
no trajectory actually touches the surface of the object, and that 
the closest approach in this figure is slightly exaggerated by the 
thickness of the line. 


the solutions are known for this specific form of third order 
polynomial 


rm = -- 


X cos 


L2 +A: 


1 -1 
- cos 
3 


3^3 


2 yi^Tfc 


27rm 


(30) 


with m = —1,0,+1, analogous to the vacuum case with 
G ^ L'^ + k. 


3.1 Circular orbits 


Circular orbits with radii occur at extrema of the poten¬ 
tial, where the first derivative vanishes. For general h this 
is 


dV^ff ^ -2L" 
dr r^ 



k 


h- —(h + i) 

T 


(31) 


Stability is determined from the sign of the second derivative 


dVeff 

dj-2 




k{h + l) 



2M 


{h 2) 


•(32) 


We summarize the behaviour of the circular paths in Table 
[T] for integer h = 0, ...,3. It is also easy to find the corre¬ 
sponding angular momentum necessary for circular motion 
in general. Setting equation m equal to zero and solving 
for the angular momentum gives 




Mr^ k 

(rc - 3M) 


'h + 1 

. rc 


h 


(33) 


This is the angular momentum that is required to keep a 
photon in a circular path around a compact object sur¬ 
rounded by a plasma distribution with a power-law density. 
The radii of the circular orbits is listed in Table[T]as Tc, along 
with the corresponding density index h. Only the h — 0 case 
admits both solutions, which are physical and become stable 


Table 1. Circular orbit radii for corresponding integer power-law 
density index h. The radial range where unstable circular orbits 
exist is also given. 


Index 

Radius of orbit 

Range 

h = 0 


3M - 6M 

h=l 


3M-4M 

h = 2 

rc = 3M 

- 

h = 3 

„ - 3M ( t 2 k ,/i-4, fe2 ,7fci.2\ 

2L^ 1 2M + Y ^ 4M^ ^ 9 M j 

SM - 3M 


when they coincide at r = 6M. For the other cases, only one 
of the roots is positive. All of the solutions have L —^ oo as 
r —^ 3M. 

Introducing terms into the effective potential through 
refractive plasma brings new behaviour to the trajecto¬ 
ries of light rays, some of which correspond analogously 
to known cases. For example, the h = 3 case is related to 
the effective potential found in the Klein-Gordon equation 
that describes the evolution of scalar field perturbations on 
the Schwarzschild space-time ( Chandrasekhar fc Detweilleij 
ll975l : [Natario fc Schiappall20oRskakala fc Visseill2010l ~). To 


proceed, we assume the energy density of the scalar field 
to be small and its influence on the space-time geometry 
negligible. Due to spherical symmetry, the radial and an- 

J ular parts o f the Klein-Gordon equation can be separated 
Zerellilll970l j. and the field written in terms of vector and 
spherical harmonics. This radial p art is the Regge-Wheeler 
equation llRegge fc Wheelerlll957lj for massless, electromag¬ 
netic (s = 1) and gravitational (s = 2) perturbations on the 
Schwarzschild background 


^ -b - I/rw) ipi = 0 


(34) 


where radial coordinates with an asterisk denote the use of 
a ‘tortoise coordinate’, 


V 


dr V r } 


(35) 


The potential governing the behaviour of the perturbation 
is 


(36) 


/ 2M\ 

+ 2M(l-s^)' 


^2 ' ^3 


with the orbital quantum nu mber I and spin given by s, 
along with the condition I > s (iLun fc Fackerelilll974l i. Writ¬ 
ing A = l{l + 1) and /3 = 1 — the solution for circular 
orbits is then 


3A4 / , 

rRw=^ A-/I + 


A" + /32 + y/SA 


(37) 


Now consider the case of a plasma surrounding a compact 
object with a density profile that drops off as 1/r®, such that 
/i = 3 in equation (1281) . The resulting effective potential is 
a direct analogy to the potential that arises in the Regge- 
Wheeler equation. The solutions for circular orbits raw and 
Tc for /i = 3 in Table [T] can be transformed into one another 
given the substitutions G —^ l{l + 1) and k — ^ 2M(1 — s^). 
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The analogy works directly for s = 0, and reproduces the 
vacuum case for s = 1 (fc = 0) as expected for photons. For 
s = 2 we must allow A: < 0 which corresponds to a refrac¬ 
tive index that has n > 1 everywhere, but is an unphysical 
density distribution for plasma. However, with this caveat 
we are able to reproduce the potential and circular orbits 
of massless scalar fields of arbitrary spin. While we work 
in the geometric optics limit, our trajectories point along 
the propagation directions of the corresponding solutions to 
the Regge-Wheeler equation. Thus, the interaction of pass¬ 
ing electromagnetic radiation with a gravitational field and 
a plasma that has a density distribution that drops off as 
1/r® is in some sense analogous to the more general inter¬ 
action of a massless boson and scalar field perturbation on 
the Schwarzschild background. 


4 EFFECTS ON PULSE PROFILES AND 

TIMING 

As a physically motivated scenario, let us consider an NS 
embedded within a plasma with bright polar emission re¬ 
gions near the surface that produce observable broad-band 
pulsations. This toy model provides an example of the signif¬ 
icance refraction can have on the self-lensing of a compact 
object. The basic model and viewing geometry used here 
were thoroughly discussed by IPechenick, Ftaclas fc CohenI 
lll983ll . and the interested reader should consult that work 
for full details. 

Let us define the observer’s coordinate system with re¬ 
spect to a plane perpendicular to the line of sight (the ‘de¬ 
tector’) which collects the rays from a distant object. The 
detector surface records the flux received from the object 
(the ‘image’) over the angles 0' and R such that an element 
of solid angle in the observer coordinate system is given by 

dfl' = sin(6l')d6i'd<(>' rs OAO'dR (38) 

under the appropriate small angle assumption for 9'. This 
angular patch is related to the impact parameter in a simple 
way by fo = D6' which gives 

dn' = (39) 

where b is given by equation (I20II and varies from 0 to 6max. 
Since the star is spherically symmetric R = cj>. Thus we can 
imagine tracing each patch of image back to some location 
on the surface of the star given by the angles A feauation ll8l) 
and (j>. In this sense A translates the two dimensional solid 
angle element of the observer to a corresponding surface area 
element on the star. To generate pulse profiles we must find 
the flux from a bright emission region at the pole of the star 
and describe the projection of this feature on the observer’s 
sky. 

The line of sight passes through the origin of the detec¬ 
tor plane and the centre of the star, taken to be the origin 
of a spherical coordinate system fixed relative to the surface 
of the star (the ‘object’ coordinates). Let a circular emis¬ 
sion region (a ‘hotspot’) be centred on the polar cap with 
an opening angle 9^, and the centre of the polar cap be de¬ 
noted by a unit vector C. Let us also define a vector along 
the line of sight, L, such that the angle between L and the 
cap centre C is called 9o. 


In general 9o is a time-dependent quantity since the 
star rotates. To express this as a function of time, we de¬ 
fine the phase 7p = fit -|- 70 where = 2tt/P with P the 
rotation period and 70 an arbitrary phase constant. Then 
the angle between the rotation axis f and line of sight is 
X = cos“^ (r ■ L) and the angle between the rotation axis 
and the centre of the hotspot is ^ = cos“^ (f • C^. The ori¬ 
entation of the cap can then be expressed with respect to 
these angles as 

9o{t) = cos“^ [cosxcos^ — sinxsin^cos7p(f)] (40) 


where the time dependence of 9o has been explicitly stated. 

In the object coordinates, an arbitrary point on the cap 
boundary is given by (0c, $), and the corresponding point 
on the observer’s detector is {9, (ph)- In general the two coor¬ 
dinate systems are rotated about the L-C plane. Using the 
transformation between the cartesian components of the two 
systems, we arrive at an expression for 4>b{9) in terms of the 
observer coordinates. 


4>h{9) = cos 


/ cos 0c — cos 00 cos 0 
V sin 00 sin 0 


) 


(41) 


for a given 0. Thus we are able to determine the boundary 
of the hotspot in a slice along 0 in the observer coordinates 
for a given orientation of the star (i.e., a given 0o) and a 
given cap size 0c. The hotspot extends over the range —4>h 
to 4>h, and we denote this range as 14(0; 0o, 0c) = 2((>b, sim¬ 
plifying the observer’s solid angle element. Determining the 
observable flux from a hotspot therefore only requires a one¬ 
dimensional integral over the impact parameter b. 

In general the definition of 4>h is more complicated than 
described above. In fact, 4>h is not defined for all arbi¬ 
trary orientations of the cap since equation m requires 
00 -b 0c < 0 < TT and 00 — 0c > 0 for example. However, the 
conditions necessary to contin uously extend cph are given in 
a vari ety of sources including Pechenick, Ftaclas fc Coheni 
(|l983l l and IPabrowski fc OsarczukI (199^ who include a 
thorough discussion for cases when multiple imaging is im¬ 
portant. To compare our results with the vacuum case we 
introduce one final substitution. In equation m we set 
X = b/M, which brings us to 


dfl =-^i4(0[a:]; 00 , 0c)a;da:. 


(42) 


To find the flux of a surface brigh t spot we must use 
the specific intensity. ISicak fc Hadrav^ (ligTSl i have shown 
that the generalization of the invariant relationship between 
specific intensity and frequency along a ray in a refractive 
medium is 


lobs 


= constant 


(43) 


where the labels obs and em denote the observed and emit¬ 
ted quantities, respectively. Taking the observation point 
to be infinitely far from the star, Wobs = ^ 00 , n-obs = 1, 
and assuming the emission point on the surface cUem = 
cUooA(i?)~^/^ and riem = n{R). By equating the observed 
and emitted intensity through this relationship we obtain 

/ 13/2 


For simplicity we take /em = rd{R)IofB,{5) where 7o is a con- 
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stant and /b a function of emission angle to accommodate 
anisotropic emission effects 

I.o,s= lofBiS). (45) 

Other choices for Jobs are also possible, such as approxi- 
mating the specific intens ity via the brightness temperature 
dLorimer fc Kramerll2005l l. The emission angle at the surface 
is a function of x since 

S = sin”^ (a;/®niax) (46) 

where Xmax = 5max/M. The flux at a given frequency oj is 

X / IofBiS[x])H{e[x]-,eo,e,)xdx. (47) 

Jo 

The result looks simil ar to that obtained by 
iPechenick, Ftaclas fc CohenI ill983l l: however, the key 
differences come from A and the range of x, since it 
depends on the maximum impact parameter which is 
now dependent on frequency. The presence of an antipo¬ 
dal hotspot simply requires a second component with 
00 —>■ TT — 00- To generate light curves we recall the time 
dependency of 0o and plot for a variety of frequency 
ratios LJe{R)/uj{R) and mass-radius relationships R/r^. The 
ratios we used start with a case that demonstrates multiple 
imaging of the stellar surface. We decrease this ratio in 
each subsequent plot to show how the pulse prohles change 
for larger ratios of radius to mass. We consider a constant 
size polar cap fixed with aperture 0c = 5° in an orthogonal 
configuration with X = C = We plot the h — 3 case 
in Fig. [4] and an analogous example with two antipodal 
hotspots in Fig.[S] We emphasize the change in morphology 
of the pulse prohles as a function of frequency by assuming 
isotropic emission, setting /b(5) = 1. The calculations are 
done in terms of the ratio u)e{R)/io{R), so the specihc value 
of k will affect only the frequency at which plasma effects 
become signihcant, and will leave the morphology of the 
calculated pulse prohles unchanged. 

We vary the angles x ^nd ^ in Fig. [6] to demonstrate 
the frequency-dependent pulse morphologies for a variety 
of viewing angles and system geometries. With t he conve¬ 
nient classihcation scheme of iBeloborodovI ll2002ll we show 
separate examples of each pulse class for antipodal conhgu- 
rations. In this scheme, class I has a primary spot which is 
always visible, and an antipodal spot which is never seen. 
In class II prohles, the primary spot is always seen and the 
antipodal spot is only visible for part of the time. Class III 
prohles are generated when the primary and antipodal spots 
are both visible for fractions of the time and class IV has 
both spots seen at any time, producing a relatively constant 
prohle. We plot the prohles for the h = 3 case which shows 
how the prohles evolve as a function of u). 

Another signihcant effect for pulse prohles is the time 
delay, which is also affected by the medium. Consider the 
arrival times of photons leaving the stellar surface and trav¬ 
elling to a distant observer 

The time delay /S.t[R) for a given impact parameter is found 


R/r =1.6 R/r =1.675 

g g 




Figure 4. Pulse profiles of a single hotspot, h = 3. The pulses are 
plotted for an orthogonal configuration with y = ^ = 7r/2, and 
frequency ratios a;e(i?)/^(i?) =0.99 (dashed Vine), LJe{R)/u}{R) = 
0.60 (dashed-dotted), LJe{R)/(j^{R} = 0.30 (dotted). The aperture 
angle is 9c = 5°. The vacuum case is plotted as a thin, solid line. 
The upper left panel has R/vg = 1.6, the upper right R/vg = 
1.675, lower left R/vg = 2.0 and lower right R/vg = 3.0. The 
pulse profiles narrow as cj(i?) approaches Ue{R). 
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Figure 5. Pulse profiles of antipodal hotspots, h = 3. The pulse 
profiles for two antipodal hotspots are plotted. Parameters are 
identical to those used in Fig. 


by comparing to the fiducial value 5 = 0. For photons emit¬ 
ted from the surface at R, we have 




Ab2\F2 


dr 


which reduces to the vacuum case when n = 1. 


( 49 ) 
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Class I 
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Figure 6. A demonstration of pulse classes for a NS with R/vg = 
3. All other parameters are as in Fig.|3 We plot plasma frequency 
ratios u}g{R)/ui{R) = 0.99 (dashed line), u)g{R)/(jj{R) = 0.60 
(dashed-dotted), uJe{R)/u}{R) = 0.30 (dotted), which is similar 
to the vacuum profile for all classes. The upper—left panel has 
X = 20°, ^ = 30°, the upper-right panel x = 30°, ^ = 60°, the 
lower-left panel has x = 60°, § = 80° and the lower-right panel 
has X = 20°, 5 = 80°. Except in class I, the profiles completely 
vanish nearest the plasma frequency. 


5 DISCUSSION 


The diverse population of NSs is defined by a wide range 
of magnetic fields, varying from the rapidly spinning mil- 
lisecond pulsars (M SPs) with fields from 10® to 10^° G 
dSacker et al.lll9^ ) to the slowly rotating magnetars which 
have s uper-critical QED scale fi elds between 10^^ and 
10^® G (iDuncan fc ThompsonllToOd ) . Generally, NSs possess 
complex magnetospheres containing p lasma that can have 
a sig ni ficant effect on p u lse pr of iles llCheng fc RudermanI 
1979: lArons &: Barnard 1 19861 : iBarnard fc Aronj 19861 : 


Wang. Wang fc Hanll2014 ) , though the exact mechanism 


sponsible for radio emission in the magnetosp here is not well 
understood (iHankins, Rankin fc Eilekll20ld ). The GJ num¬ 
ber density provides a lower limit for the charge density in 
the neighbourhood around an NS, and is given by 


Ngj = — -COS 4 ; 


7 X 10 


— 2 Bz 


(50) 


27rec P 

in units of particles per cm®, with P the rotation period in 
seconds and Bz (G) the field co mponent along the rotatio n 
axis in the approximation above (iGrrldreich fc Julianll 196911 . 
Since Nqj is a lower limit we use the proportion ality k to 
discu ss the magnetospheric density No = kNqj llSturrockI 
11973). This expression gives a physical motivation for 
the /i = 3 plasma density index discussed throughout 
this work. The ray-tracing discussed in Section [5] can be 
used for rays passing through a pulsar wind, introducing 
relativistic and plasma effects into models of binary eclipses 
for a double NS or NS-black hole binary that is nearly 
edge on with respect to the observer. In the case of an 
aligned rotator we expect a wind density proportional to 
1 /r, such that h = 


]__(l22^22SUl2Si=fl=pl2s_^ 


]j99^ _ iHar ding. Gontopoulos fc Kazan^ _ 

Kirk. Lyubarsky" fc PetTi l2009l : Tong. Xu. Song fc 


idtl 


19991: 


Qiad 


l2013h . shown in Fig. [2] It is plausible that both lensing 
and plasma may affect the eclipse mechanism for as-yet 
undiscovere d compact binary systems -with a n edge on 
orientation llKeane et al. li^; IShao et al.| [20li). 

Normal radio puls ars show a rich varie ty of frequency- 
dependent behaviours llChen fc Wang|[20l4l . A proportion¬ 
ality constant 10 ® < k < 10 ® is needed to describe the dense 
environment surrounding these objects (Ilstominll200ll) . Gen¬ 
erally the simplification of cold plasma is unrealistic since 
matter within the light cylinder is strongly magnetized in 
corotation with the NS and highly relativistic. Particles 
moving along the field lines have a Lorentz factor 7 = 
EfmgR, giving an effective plasma frequency 


tUeff 


OJe 


(51) 


where the right-hand side is summed and averaged 
over the positron and ele ct ron d ist ribution functions 
([Ge dalin. Melrose fc GrumarJ 1 19981 : iMelrose fc GedalinI 
Il999lb The relativistic case giv es a much lower cutof f 
frequenc y, we ll below 100 MHz ll Melrose fc GedalinI Il999l : 
iMelrosel 20001) . Thus, despite the low frequencies involved 
the interplay between GR and plasma eff ects in pulsars may 
potentially be within ob s ervable ranges ( Pgggvet_alJ I 2 OO 6 I : 


iMalov fc Malofe"^ l2010l : IStappers et al. 2011^ . It is also 


significant to note the existence of a subset of radio pulsars 
with pulse widths th at broaden with increasing frequency 
lichen fc Wand[2014l) . also seen in Figs. [4] and [5] for the toy 
light curves generated here. 

Due to their fast rotation rate, the radio emission re- 
gion is thought to be closest to the NS surface in MSP s 
lAlpar et al.l Il9^ : iBhattacharva fc van den Heuveli Il99ll ) . 
If a low-frequency component of the radio emission from 
M SPs does originate sufficiently near th e surface (ie, r/R < 
10: lGil fc Kiialjll993l:liaak fc Gillll998l) . it may be plausibly 
affected by both relativistic and plasma effects. The signifi¬ 
cance of rotation for a quickly spinning MSP, especially due 
to induced oblateness of the surface, becomes increasingly 
important for the pulse profiles and time d elays, making 
the S chwarzschild case a first approximation llCadeau et al.l 

I2OO7I) . 


At the opposite extreme from MSPs, anomalous X-ray 
pulsars (AXPs) and soft gamma repeaters (SGRs) appear 
to be NSs with magneti c fields above the QED critical field 
Bqed = 4.4 X 10^® G llMereghettil l201 3l). most commonly 
descr i bed by the mag n etar model llDuncan fc Thompson! 
I 1992 I : IPaczvns^ 1 1992 : lUsovI Il992l ). The lack of persis- 
tent radio emission abov e 500 MHz lllstomin fc SobvaninI 
I2OO7I : liTazarus et al.ll20l3 ) and the transient, highly-variable 
radio emission observed from a subset of these objects 
llGamilo et al.l l2007l : iMalovI l2014l ) suggests that the den¬ 
sity of plasma in the magnetospher es of magnetars is 
lower than in norm al radio pulsars lllstomin fc SobvaninI 
I2OO7I : iMalovI I2OI4I ). While gravitational effects are ex- 
pecte d to dominate photon trajectori es outside of magne¬ 
tars llHevl fc Shavivl]2002l : iKimi [ioi^ . it is generally be¬ 
lieved that plasma effects should also become im portant 
deep within magnetar atmospheres llLai fc Hdl2003h . A par¬ 
ticularly interesting alternative to the magnetar model de- 
scribes AXP and SGR behaviour in terms of a quark star 
llOuved et al.l l2004l : iLeahv fc OuvedI l2009fl . The environ¬ 
ments around these objects are expected to be complex, 
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and some models explain magnetar-like behaviour through 
accr etion from fail-back material on to the surface of the ob¬ 
ject (jOuved. Leahy fc Niebergalll2007t l. Interest in the quark 
nova scenario has grown considerably due in part to the 
NuST AR observations of the C assiopeia A supernova rem¬ 
nant llGrefenstette et ^ l2014h which detected a discrep¬ 
ancy in the distribution of titanium-44 and iron-56 abun- 
dances that can be explained by the quark nova scenario 
dOuved, Leahy fc Koninal2015l ). Considerations of more re¬ 
alistic dispersion relations coupled with GR may be relevant 
in the dynamic environment around an energetic, highly 
compact quark star. 


6 CONCLUSIONS 

Ray-tracing through the Schwarzschild space-time includ¬ 
ing plasma demonstrates light ray paths where refraction 
can dominate over gravity, scattering incoming rays. Figs. 
[2] and [3] show the effects for density power-law index h — 1 
and h — 3, respectively. We have explored the orbits dis¬ 
cussed in Section [3] with both the integral formalism (equa¬ 
tion mi and by solving the system of differential equations 
given by equations [TT] and [T^] directly using a fourth- order 
Runge-Kutta method. The integral approach is valuable for 
analytical results like the pulse profiles given in Section U 
while the RK4 scheme is suited to handling general cases 
when spherical symmetry is not assumed (for example, with 
the Kerr metric), and will be useful for more complicated 
plasma distributions. 

In Section [331 we found the radii of circular orbits for 
arbitrary h and derived the turning points for trajectories 
that do not reach an observer at infinity. Depending on the 
details of the density distribution the plas ma frequency can 
act a s an effective photon mass (h = 0; iKulsrud fc LoebI 
I 1992 II . an additional contribution to the angular momentum 
{h = 2), or reproduce the potential of a massless scalar 
field perturbation on the Schwarzschild background with a 
judicious choice of scaling constant {h = 3). 

The interaction between a refractive medium and 
gravitational lensing has b een considered at a variety 
of scales in the literature ([Broderick fc BlandfordI I 2 OO 3 I : 
iBisnovatvi-Kogan fc Tsupkcill2009l : lEr fc Madl2014li . The ef¬ 
fect of plasma on the formation of gravitational lensed 
images of background sources will be difficult to de¬ 
tect due to the small angular size of the images. Ob¬ 
serving frequency-dependent magnification from transient 
microlensing events requires mu l ti-wav eband monitoring 
llMorozova, Ahmedov fc Tursuno^ l2013l i . In this work we 
investigated another scenario in which dispersion and gen¬ 
eral relativistic effects may overlap: the pulse profiles gen¬ 
erated by a compact object embedded in a continuous and 
smoothly varying refractive medium. We made use of the 
Schwarzschild metric to study the frequency-dependent ef¬ 
fects, but a similar anal ysis can be done using the Reiss- 
ner N ordstrom and Kerr llMorozova. Ahmedov fc TursunovI 
I 2 OI 3 II solutions. Both charge and rotation influence the 
space-time ge ometry and therefore wou ld produce unique 
lensing effects llChakrabortv fc Senll2014 f. This work shows 
that inclusion of even the simplest plasma effects generates 
achromatic behaviour in pulse profiles and novel trajectories 
for light rays near massive compact objects. 
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